NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



( It C ' J 


ASRL TR 196-3 

(NASA-Cn-1to6233) TESTING AND EVALUATION OF 881-29135 

A STALL” F’-UT1 ER-SUPPRESSXON SISTEfl FOB 
ftELICOPTBR BOTOBS USING 

iNDI VIDUAL-BLADE-CX'HTBOL (Massachusetts Uacias 

lust, or Tech.) 62 p HC A05/MF A01 CSCL 01C G3/U 8 32 902 

TESTING AND EVALUATION OF A STALL-FLUTTER-SUPPRESSION 

SYSTEM FOR HELICOPTER ROTORS USING 

INDIVIDUAL-BLADE-CONTROL 



Distribution of this report is provided in the interest of 
information exchange. Responsibility for the contents re- 
sides in the author or organization that prepared it. 


Prepared under research Grant No. NSG-2266 by 

Aeroelastic and Structures Research Laboratory 
Department of Aeronautics and Astronautics 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 

for 

AMES RESEARCH CENTER 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
MOFFETT FIELD, CALIFORNIA 94035 


I 


ii 

TESTING AND EVALUATION OF A STALL FLUTTER SUPPRESSION 
SYSTEM FOR HELICOPTER ROTORS USING 
INDIVIDUAL BLADE CONTROL 

by 

TODD RANDALL QUACKENBUSH 


ABSTRACT 


The development and testing of a feedback system designed 
to alleviate the violent blade first torsion mode oscillations 
associated with stall flutter are described. The system, based 
on previously developed M.I.T. Individual-Blade-Control hard- 
ware, employs blade-mounted accelerometers to sense torsional 
oscillations and feeds back rate information to increase the 
damping of the first torsion mode. A linear model of the blade 
and control system dynamics is developed and is used to give 
qualitative and quantitative guidance in the design process as 
well as to aid in analysis of experimental results. System 
performance in wind tunnel tests, both in hover and forward 
flight, is described, and evidence is given of the system's 
ability to provide substantial additional damping to stall- 
induced blade oscillations. 
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SECTION 1 


INTRODUCTION 

As increased demands on helicopter performance have pushed 
machines to higher values of blade loading and advance ratio, 
one persistent problem for the designer has been the transient 
aeroelastic instability known as stall flutter. This phenomenon 
has been extensively studied in a variety of other works (Refs. 
1-5, to name just a few) and a comprehensive discussion of its 
sources and effects is not necessary for present purposes. 
However, a brief summary of the salient points is helpful for 
posing the design problem dealt with herein. 

It has been well documented that an airfoil oscillating 
rapidly in pitch is able to operate transiently at angles of 
attack considerably in excess of its static stall angle without 
flow separation taking place. However, at sufficiently high 
angles of attack, the airfoil stalls, though this so-called 
dynamic stall differs considerably from conventional static 
airfoil stall. As shown by Ham, Ref. 5, dynamic stall is 
characterized by the loss of leading edge suction and the 
subsequent movement of a large negative pressure disturbance 
aft from the leading edge, a movement which generates strong 
nose-down pitching moments on the airfoil, with proper 
combinations of airfoil mean angle of attack, amplitude of 
motion, and reduced frequency, this stalling phenomenon can 



cause an aerodynamic moment hysteresis which can lead to a net 
influx of energy to the airfoil's pitching motion. 

The application of these findings for airfoil dynamic 
stall to helicopter rotor blades is relatively straightforward. 
Figure 1 shows a typical azimuthal angle of attack distribution 
for a rotor blade in forward flight. The high angles of attack 
on the retreating side and the rapid pitch angle variations 
caused by cyclic pitch inputs strongly suggest that rotor 
blades under certain conditions should be susceptible to the 
same stall-induced oscillations observed in airfoils. Both 
experiments and flight experience has shown that this is indeed 
the case; for certain combinations of blade torsional natural 
frequency, blade loading, and advance ratio, the spanwise 
integrated effect of dynamic stall is to feed energy into blade 
torsional motion, particularly the first torsion mode. This 
motion is generally only transiently unstable and damps out 
rapidly as the blade swings around the azimuth toward the 
advancing side. However, even the one or two cycles of blade 
motion that do occur are sufficient to put extreme loads on 
the rotor control system (see Fig. 2); the fatigue life of 
rotor pitch links can thereby be considerably reduced. 

Within the restrictions of conventional swash-plate control 
systems, methods to alleviate this problem all have considerable 
drawbacks. Increasing solidity would reduce blade loading and, 
thus, susceptibility to stall, but this would penalize overall 
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helicopter performance. Restricting the flight envelope of the 
machine to low advance ratios is undesirable for the same reason. 
Use of airfoils with more benign dynamic stall characteristics 
is possible , but this further complicates the already very 
involved rotor airfoil design process. 

Applying Individual-Blade-Control (IBC) techniques to 
this problem offers a possible solution. Reference 6 showed 
that appropriate feedbacks to a position control servo govern- 
ing blade pitch motion could help reduce undesirable blade 
motions due to low-frequency gust inputs. It was felt that 
similar methods could be applied to alleviate the violent 
torsional motions associated with stall flutter. To understand 
the overall concept that was employed , consider again for a 
moment the mechanism which drives the stall flutter oscilla- 
tions. As noted previously, at high blade angles of attack 
and certain reduced frequencies, aerodynamic moment hysteresis 
causes a net input of energy to blade torsional motion, so that 
any small blade oscillation grows with time. Such a situation 
is reminiscent of simple oscillating systems operating with 
negative damping. Indeed, even though stall flutter of rotor 
blades is in reality a result of aerodynamic forcing, it can 
be conceived of as a phenomenon caused by a once per revolution 
variation in the effective damping of the blade in pitch. On 
the advancing side, the blade experiences strong positive 
damping at low angles of attack, but on the retreating side 
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the affective damping can temporarily become negative* leading 
to the oscillations described above. Figure 3 gives an idea 
of the variation of the effective damping function with blade 
operating condition. 

An effective stall flutter suppression system* then* would 
be one which could eliminate this one-per-rev excursion into 
negative damping. One way to achieve this end which is 
suggested by classical control theory is to provide a pitch 
rate feedback from the blade to the pitch control servo. The 
details of the rationale for this concept* its implementation* 
and the results of experiments based on it are given in the 
following sections. 
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SECTION 2 

MODEL DESIGN AND I.B.C. HARDWARE 

2.1 Overall Modal Construction 

The model used here to test the proposed stall flutter 
suppression system was identical in most particulars to that 
used in Ref. 6. A D.C. servomotor serving, through a series 
of linkages, as a blade pitch position control system was 
mounted on the rotor shaft. The test rotor used only a single 
blade, with a NACA 0012 section, 21.2 inch span, and a two 
inch chord; further details on the blade are given in Table 1. 
The blade was attached to the rotor hub by means of an aluminum 
fork which in turn was connected to spherical bearing mounted 
in a ball and socket assembly; thus, the blade's flapping, 
lagging, and feathering motions all took place about the same 
point. A steel flexure instrumented with strain gauges was 
attached to the blade to sense pitch angle. 

Two "dummy blades" in the form of lengths of threaded 
5/8" steel rod were also attached to the rotor hub. Each rod 
had adjustable counterweights which were used to achieve dynamic 
balancing during rotor operatio* Two symmetrically mounted 
counterweights were also attached to the shaft to balance the 
mass moment of inertia contribution of the active motor. 

Photographs at the blade and control system hardware are 
shown in Figs. 4-6. Further details of the construction of the 


actuation system are given in Ref. 6 and will not be repeated 
here. 


2.2 Introduction of Blade Torsional Flexibility 

Since the primary aim of this experiment was to design 
a system to control the first torsion mode of the rotor, it 
was necessary to ensure that the frequency of the mode was 
within the bandwidth of the servomotor. For full scale rotor 
blades, Wg is usually of the order of 5fi to 7fl, or about 30-40 
Hz for most helicopters. Unfortunately, values of Wg for small, 
relatively stiff model blades such as the one employed here 
are invariably much higher than for full scale blades; in this 
case, it would have been very difficult, because of the high 
value of oj Q in the test blade, to induce the blade itself to 
flutter. Even had that been achieved, due to restrictions on 
the servo bandwidth (approximately 40 Hz) , control of these 
oscillations would have been difficult. 

Given this, it was necessary to introduce torsional 
"softness" into the control system artificially. Two different 
methods to achieve this end were used here: 

First, adjustment of rate and position feedbacks in the 
servo position control system were made to produce an "electronic 
spring" at the blade root. For the actuation system described 
in Ref. 6, the connection between the servomotor and the blade 
is essentially completely rigid. However, note that for this 
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case the transfer function of the position control system can 
be written as follows (neglecting pitch flexure high frequency 
dynamics) : 


_ zrsi l.O 

V ~~ s' + (.iSl KiV * WZOO H § 


(See block diagram 
in Fig. 7) 


( 1 ) 


As is evident from the above , if were reduced to very 
small values, the actuator would closely resemble a very lightly 
damped torsional spring, with the spring’s natural frequency 
determined by K Q . Such a situation is a reasonable simulation 
of an actual rotor blade, since, for full-scale blades, most of 
the torsional "softness" originates in the control system itself. 

Second, additional "softness" was introduced by inserting 
two leaves of spring steel (dimensions 3-1/2" x 1-1/8" x .02") 
between the blade mounting fork and blade itself (see Fig. 8) . 

As shown in Fig. 8, the leaves were installed so that they lay 
parallel to the plane at the blade when at rest. This modifica- 
tion (which added 3.5" to the blade radius) was achieved by 
bolting one end of the leaves to the mounting fork which 
originally clamped directly onto the blade (with the aid of 
steel filler plates in the fork) ; the other ends of the leaves 
were secured to the blade with the aid of two 2-1/4" x 1-1/4" 
x 3/16" aluminum clamping plates built for this purpose (steel 


fillers were again used between the clamping plates) . Details 
of the "tuning" of this arrangement to give the desired are 
given in Section 4. 

2.3 Pitch Rate Signal Extraction 

As noted in Section 1 , the principal aim of this experi- 
ment was to increase the damping of blade torsional oscillations 
by supplying a pitch rate feedback to the position control 
system. It was decided to extract the pitch rate signal by 
first obtaining a pure pitch acceleration signal from acceler- 
ometers mounted on the blade and then feeding this signal to an 
integrator. 

Appendix 1 discusses in detail the signals sensed by 
accelerometers mounted a distance a away from the pitching axis 
of the blade. As noted in that Appendix, any single acceler- 
ometer mounted away from the pitching axis will sense a component 
of centrifugal force (often called "propeller moment") which is 
proportional to pitch angle. This difficulty can be overcome 
if the signals from two separate accelerometers mounted of an 
equal distance from the pitching axis but oriented with a 90° 
separation between them are summed; the result is a signal 
purely proportional to 0 . 

Unfortunately, the above result is strictly true only 
if the flapping and lagging degrees of freedom are neglected. 
Unless the accelerometers are placed such that they lie in the 
plane formed by the lagging and flapping axes, components of 
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centrifugal force proportional to lag angle £ and flap angle 3 
will enter into the signal sensed by the accelerometers. This 
situation has the potential for causing difficulties with the 
integrated feedback signal, since an ideal integrator would 
apply an infinite d.c. gain to any steady-state components in 
the £ and 3 signals. Fortunately, the actual integrator used 
here was non-ideal, having a roll-off at low frequencies (below 
approximately 0.5 Hz) (see Fig. 10) which would eliminate any 
steady-state signals. However, blade lag motion characteristi- 
cally also involves a low frequency component at approximately 
0.2Q to 0.40; to minimize the possibility that the integrator 
would amplify these low frequency signals (at approximately 1.3 
to 2.6 Hz for the rotor speeds used here) and distort the 
feedback, the brackets holding the accelerometers were built 
such that they were rigidly attached to the blade but were 
located in the plane of the lagging and flapping axes (see 
Fig. 9) . These brackets were sized to be used with the 
"electronic spring" blade configuration i.e., when the spring 
steel inserts were not present between the actuator fork and 
the blade. When these inserts were used, it became impractic- 
able to have the accelerometers lie in the plane of the flapping 
and lagging axes; the same accelerometer brackets shown in 
Fig. 8 were used in this case without apparent adverse effect 
on system performance. 
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The accelerometer installation shown in Fig. 3 will 
inevitably sense a strong 1ft signal if a pilot pitch input 
(cyclic) is used. This circumstance will not interfere with 
the intended purpose of this feedback control system (i.e., 
providing additional damping to transient disturbances in 
torsion) since the 1ft is a steady input; a rate feedback on 
such a signal will merely introduce a phase lag which can be 
compensated for elsewhere in the pitch control system. 
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SECTION 3 

THEORETICAL ANALYSIS OP CONTROL SYSTEM 

To see the fundamental concept behind the flutter 
alleviation system studied here, first note the system block 
diagram in Fig. 7. This diagram illustrates the "electronic 
spring" case, in which torsional softness is obtained by 
reducing Kg to low values and then adjusting Kg to obtain the 
desired "spring" undamped natural frequency. The open loop 
transfer function of this system for ft * 0 is given by 



P x = -182.7 + 416. Oj p 2 * -164.2 + 328. 5j 


The above poles p^ and p 2 are obtained from Fig. 7 for 
Kg * .234 and Kg = .05, with the pitch flexure dynamics included. 
As the root . locus diagram (Fig. 11) shows, the inclusion of the 
pitch flexure dynamics (discussed in Appendix 2) considerably 
changes the simple system model one would have expected were 
these dynamics neglected, as they were in Eq. 1. However, 

Fig. 11 also shows that the feedback nevertheless achieves the 
fundamental aim, i.e., increasing the damping of the blade 
"first torsion mode". It is clear, though, that the range of 
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permissible K p values is limited, since large K p drives the 
flexure pole unstable. 

For the case in which flexible leaves are inserted between 
the actuator fork and the blade, the dynamics of the system 
become somewhat more involved, since in this case the blade 
pitch angle can differ from the motor shaft command pitch angle, 
unlike the earlier case in which the shaft/blade connection was 
essentially rigid. The changes effected in the system equations 
of motion and in the overall system dynamics by this additional 
degree of freedom are shown in detail in Appendix 2. For this 
case, the effect of including the pitch flexure dynamics was 
negligible. The open loop transfer function for the "mechanical 
spring" system for ft * 0 is then 

( W 1 ) (£».<> *•) (i- ’4)(i-%)( 14 KJ (3 > 

f,= -i.-ir *111.1 j 

The block diagram of the system is shown in Fig. 12 and 
the root locus diagram is shown in Fig. 13 for = .528 and 
K* = .19. These values of K Q and were chosen since they 
gave what appeared to be an acceptably large range of values in 
which the rate feedback could succeed in keeping the blade 
oscillations stable; differen values of feedback gain, for 
example Kq * .528, = .08, would produce a root locus 
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diagram such as in Fig. 14, in which any but very small values 

of K_ drive the blade oscillations unstable. 

F 

Note that the system dynamics detailed above neglect the 
effects of aerodynamic damping (be it positive or negative) on 
the blade. Inserting a term D^6 in the equations of motion of 
Appendix 2 to account crudely for aerodynamic damping shows that 
the effect of aerodynamics is to move the blade oscillation pole 
of Fig. 13 to the right with negative D^, while oscillation 
frequency is nearly constant; for sufficiently large, negative 
values of D A the pole is driven into the right half-plane, but 
only transiently, since and, thus, pole location vary with 
azimuth. For proper choices of K» and K p though, blade oscilla- 
tions can be stable all around the azimuth (again see Fig. 13) ; 
nor should the value of K p required to stabilize the pitch 
oscillations be excessive, since the effective "negative damp- 
ing ratios" associated with stall-induced instabilities are 
rather small (C ef £ = -0.1 at a maximum) . 

Figure 13 illustrates two interesting aspects of the 
model developed here. First, the root locus diagram indicates 
that the blade oscillation frequency will increase substantially 
as K p is increased. Second, the analysis predicts that the 
oscillation will be stabilized for only a certain range of K p 
values and will develop a relatively high-frequency, potentially 
unstable oscillation if K p is increased sufficiently. As will 
be discussed in Section 4, both qualitative and quantitative 
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agreement was found between the predictions made above (and in 
Appendix 2) and actual experimental results. 

Finally, as is evident from the above discussion, little 
in the way of detailed aerodynamic or structural analysis was 
done preparatory to the experiments described herein; this 
circumstance came about fundamentally because it was felt that 
such analysis was impracticable and unnecessary given the 
intended scope of the present work. A number of detailed and 
complex analyses of stall flatter and its effects on rotor 
blades have already been performed with powerful computational 
tools and elaborate models (see, for example. Refs. 7,8). It 
was felt that since the system to be designed and tested here 
was to be a first-pass, proof-of-concept effort, the simple 
dynamic and aerodynamic models used above and in Appendix 2 
were adequate to the task. This evaluation was borne out by 
the results presented in Section 4. 
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SECTION 4 

WIND TUNNEL TESTING 

4 . 1 Experimental Set-Up 

Testing of the I.B.C. stall flutter suppression system 
was performed in the M.I.T. Wright Brothers Wind Tunnel. The 
7' x 10' test section contained two vertical trunnions which 
supported the rotor shaft in a horizontal attitude. This 
orientation, which caused the rotor to rotate in a vertical 
plane, was a result of the mounting requirements of the previous 
series of I.B.C. gust alleviation tests (Ref. 6). One conse- 
quence of this orientation was to introduce a one-per-rev 
gravity pulse into the accelerometers used in the control system; 
however, the magnitude of the pulses was sufficiently small 
that no adverse effect on system performance was expected or 
observed . 

The rotor was driven by an external hydraulic motor. The 
shaft was equipped with slip rings to provide power to the 
servomotor and to extract data from the various sensing elements. 
On-line data extraction was accomplished using software pre- 
viously developed by other members of the I.B.C. project team. 

4.2 "Electronic Spring" Configuration Tests 

The first series of tests involved four experiments using 
the "electronic spring" configuration. For those tests, as 
noted in Section 2, a rigid motor shaft-blade linkage was used, 
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and torsional softness was introduced by adjusting rate and 
position feedbacks in the servo (see Fig. 7) . For the values 
chosen (K Q * .234, K£ * .05), the effective damped natural 
frequency of the blade in pitch was 53 Hz. 

First, two static tests were conducted with the blade 
hanging vertically downward. A pulsed voltage signal was fed 
into the pitch control system to excite blade motion. The 
signal consisted of a series of 4.0 ms pulses at 6.0 Hz with 
an amplitude of 0.75 volts (equivalent to 3° of pitch). Blade 
pitch rate and acceleration response to this input for the open 
loop case are shown in Fig. 15. Next, the control loop was 
closed by applying the pitch rate feedback to the servomotor 
input; the potentiometer setting in the feedback loop was 
0.54. For the same excitation as the previous case, the pitch 
rate and acceleration traces shown in Fig. 16 were obtained. 
Comparing the time traces in Figs. 15 and 16, one can see that 
the damping of the rate response for the closed loop case is 
higher than that in the open loop case; the effective damping 
ratio for the latter is 0.13, as opposed to 0.23 for the former. 
Note that the closed loop case repsonse contains a substantial 
high-frequency component that is not evident in the open loop 
trace; this is the pitch flexure mode, which moves to lower 
damping and higher frequency as K p is increased (see Fig. 11). 
The amplification at this high frequency component is not 
desirable nor, though, is it material to the current aim of 
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adding damping to the blade "first torsion" mode. As Fig. 17 
(which gives the Fast Fourier Transforms of the rate responses 
for the open and closed loop cases) shows , blade response to 
the pulsed input has been substantially decreased in the vicinity 
of the original "first torsion" frequency of 53 Hz. 

A similar experiment was then performed with the same K g. 

Kg and input excitation , nowever , in this case the rotor was 
spun at 6.2 Hz; the aim here was to ensure that the control 
system would operate properly with a non-zero propeller moment 
being sensed by the accelerometers. Spinning the rotor had no 
apparent adverse effects on system performance in this configura- 
tion , and a comparison of rate and acceleration traces as well 
as the Fourier-transformed rate responses for the open and 
closed loop cases is shown in Figs. 18, 19 and 20. Here, the 
closed-loop pitch rate response has an effective damping ratio 
of 0.25, as opposed to 0.18 for the open loop case. Note also 
that the effective damped natural frequency of the rotating 
open loop case appears, from Fig. 20, to be approximately 47 Hz, 
as opposed to 53 Hz for the non- rotating case. This is not 
surprising, since rotor rotation causes aerodynamic deunping to 
be added to the blade/servo system, driving the blade pole of 
Fig. 11 farther into the left half-plane, which leads to a 
lower damped natural frequency, given the pole's trajectory. 
Despite this difference, substantial reductions in blade rate 
response around the "first torsion" frequency were observed, 
just as in the non-rotating case. 
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Though these experiments were not realistic in the sense 
that they relied on electronic rather than aerodynamic excita- 
tion, and in the sense that the dynamics of the "electronic 
spring" used here only vaguely resembled those of a real rotor 
blade, nevertheless they demonstrated that the proposed control 
system had the ability to increase the damping of high frequency 
(approx. 50 Hz) oscillations in pitch (equivalent to oscilla- 
tions in torsion for this simple dynamic system). Furthermore, 
these results had considerable qualitative and quantitative 
agreement with the theoretical predictions of Appendix 2 and 
thus gave some grounds for confidence in the analytical methods 
used. For example, for the non-rotating case. Appendix 2 
predicts damped natural frequencies of 52 Hz and 66 Hz for the 
blade and flexure oscillation modes, respectively, for the open 
loop case; the values obtained from experiment were 53 and 71 
Hz , respectively . 

Unfortunately, it was also found that the blade in the 
"electronic spring" configuration was not easily susceptible 
to stall flutter, despite vigorous attempts to excite the blade 
both in hover and in forward flight. The determination was 
therefore made to switch to the mechanical spring configuration 
described in Section 2. Previous tests at the M.I.T. VTOL 
Technology Lab had shown that this configuration was susceptible 
to stall flutter. 
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4.3 Teats with Mechanical Torsional Spring 

After the installation of the mechanical spring described 
previously# a series of tests were run with the blade static# 
in hover# and in forward flight at two separate advance ratios. 
The static tests were designed to reestablish the fundamental 
workability of the proposed control system given the substantial 
change in system dynamics introduced by the mechanical spring. 
Before these tests were conducted# a choice had to be made 
concerning the number of steel leaves to be inserted between 
the actuator and the blade# since the number of leaves deter- 
mines natural frequency of the torsional spring. One leaf 
(installed as shown in Fig. 9) yielded an ui Q of approximately 
17 Hz# while two leaves yielded 31 Hz# and three leaves 42 Hz. 
The two- leaf configuration was chosen since the frequency fell 
within the bandwidth of the servo and also provided sufficient 
flapwise stiffness to avoid interference between the blade and 
the shaft support trunnions during rotor run-up. 

The static tests proceeded much as did those with the 
electronic spring configuration#* the blade was hung vertically 
downward from the shaft while an external pulse train of 
frequency 6.7 Hz was introduced to excite blade motion. For 
these tests# and for all succeeding tests# three separate 
levels of Kp (0.0# 0.13# and 0.26) were established as bench- 
marks which covered the significant range of system performance 
and which would be used as operating points in those cases 
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where circumstances permitted. The rate and acceleration 

responses of the blade to the identical pulse input for these 

three levels of K p are shown in Figs. 21 , 22 , and 23. Graphic 

representations of the Fast Fourier Transform of the pitch rate 

signals for K - 0.0 and K * .26 are shown in Fig. 24. Both 
r P 

of these sets of results reflect the trends predicted in 

Section 3; the damping of the blade oscillations of the blade 

"first torsion" node is increased with increasing K p , while 

the frequency of the oscillation increases. Note that the 

values of K p used here correspond to the indicated points on 

the root locus diagram of Fig. 13. 

The values of K are multiplied by the appropriate 

constants (see Fig. 12) to yield the total feedback gain K p . 

The overall open- loop gain which is used in generating 

Fig. 13 is obtained by multiplying K p by the forward loop gain 

of .127 (see Appendix 2). The values of K_ corresponding to 

K p = .13 and K p = .26 are, respectively, .000135 and .00027. 

As seen from Fig. 13, these values of correspond to blade 

oscillations with c = .06 and damped natural frequency of 208 

rad/sec (for K p = .13) and £ * .085 and damped natural frequency 

of 224 rad/sec (for K p = .26). These values can be compared to 

values of .15 and 212 rad/sec, respectively, for K ■ .13, and 

P 

.24 and 233 rad/sec, respectively for K p = .26, that were 
observed experimentally. These results suggest that there was 
considerably more damping present in the system than the model 
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of Appendix 2 assumed. This is not too surprising, since only 
the mechanical friction of the motor itself was included in 
the model; many other possible sources of friction (e.g. gear 
meshing, linkage friction) doubtless existed but were difficult 
to include in a linear model and so were neglected. The 
relatively close agreement of the frequency predictions, 
though, was encouraging, as was the fact that the model predicts 
that the system will go unstable at high gain at 70 Hz, which 
agreed reasonably closely with the 65 Hz observed in experi- 
ments. 

Unfortunately, the usefulness of the model in Appendix 2 
is limited to the static rotor cases. The time-varying 
aerodynamic damping in the hover and forward flight cases 
introduces sufficient additional complications that detailed 
predictions with this model, which neglects aerodynamic effects, 
become invalid. However, as will become apparent in the follow- 
ing discussion, the overall functioning of the system is not 
impaired by this circumstance. 

To summarize , then, the effective damping ratio of the 

K =0.0 case was .02, while for K * .13 it was .15, and for 
P P 

K = .26 it was .24. From Fig. 24, one also notes that the 
P 

pitch rate component at was decreased dramatically relative 

to the K =0.0 case with K = .26. These results again show 
P P 

the strong potential of this system for increasing the damping 


of torsional oscillations in the 30-40 Hz frequency range. 
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The next step in testing was to run three cases with the 
rotor operating in hover and with external pitch excitation 
supplied. Ordinarily, it is possible (see Ref. 5) to induce 
stall flutter in hover by raising collective pitch until the 
blade stalls and provides the necessary excitation. However, 
the pitch linkage on the model rotor had not been specifically 
designed to operate at high collective, and the maximum settings 
available (16 0 -17°) were insufficient for the deep blade stall 
needed, so once again resort was made to external excitation. 

For these hover tests, the rotor was operated at 15° 
collective in order that aerodynamic damping would be low. (see 
Fig. 3). Rotor rotation frequency was 6.7 Hz, and the frequency 
of external excitation was either 6.7 Hz or 3.4 Hz, while the 
magnitude of each input pulse varied from case to case (though 
it was never more than approximately 2° or 0.63 volts). A 
comparison of the pitch rate responses for K * 0.0 and for 

r 

K =0.26 are shown in Fig. 25. (Note that excitation frequency 
P 

was only 3.4 Hz for K = 0.0; given the low aerodynamic damping, 

P 

any more strenuous excitation might have generated a divergent 
oscillation.) The effective damping ratio for the K p = 0.0 is 
0.03, while for K p = .26 it is 0.14. A comparison of the FFT 
breakdowns of the pitch acceleration response for the two cases 
in Figs. 26 and 27 shows once again the large decreases in 
acceleration components in the vicinity of = 31 Hz. Also 
evident when K p is increased, though, is an amplification of 
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pitch acceleration response at frequencies higher than (at 
approximately 40 Hz or about 6(2) . One possibility is that this 
shift is a result of the "stiffening" effect predicted in Section 
3 and noted in the description of the non-rotating tests above 
i.e. that increases in IT tend to cause w a itself to increase. 
However, such inferences are complicated by the fact that signifi- 
cant aerodynamic damping is present in the rotating tests and, as 
noted previously, such effects can also alter the blade oscilla- 
tion frequency. Hence, in the absence of a detailed aerodynamic 
model, it is difficult to isolate the cause of this frequency 
shift. Nevertheless, the overall aim of demonstrating that high 
frequency oscillations could be suppressed was achieved. 

To complete the testing series, two sets of forward flight 
cases at moderate advance ratio were run. In these cases it was 
found that stall flutter could be brought about and, thus, resort 
did not have to be made to external pitch excitation. For the 
first set of runs, rotation frequency was 6.7 Hz, tunnel speed 
was 20 mph, and hence advance ratio was 0.30. As shown in the 
top of Fig. 28, root collective pitch was 10° and a one-per-rev 
cyclic pitch signal of +9.0° was superimposed, with the maximum 
pitch angle reached at ip - 270°. The lower trace of Fig. 28 shows 
the pitch rate response of the blade; superimposed on the strong 
one per rev is the higher frequency oscillation caused by dynamic 
stall, which shows up as the slight "lumpiness" in the pitch rate 
curve. The effect is also evident in the pitch acceleration 
response in Fig. 30. 
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The effect of increasing from zero is shown in Figs. 

29 , 31# and 32. With progressively higher values of K p , the 
" lumpiness" of the pitch rate curves smooths out# representing 
a diminuation of the stall-induced oscillations. The acceleraion 
response reflects the same trends# with the high frequency oscilla- 
tions decreasing in magnitude? the FFTs of the pitch acceleration 
response shows that = .13 causes a reduction of approximately 
33% in the component at uj 0 over the K p * 0*0 case, while K p * .26 
leads to a reduction of approximately 67%. 

A second set of forward flight experiments was performed 
at essentially the same collective, cyclic, and wind speed 
settings, though the advance ratio was increased to 0.33 by 
decreasing to 6.1 Hz. This flight condition induced a substan- 
tially stronger stall flutter oscillation than in the previous 
case, so much so that, as seen in Fig. 33, K could not be reduced 

r 

to zero, only to .02, lest a violent oscillation develop. Even 
with this stronger excitation, though, settings of K * .13 and 
.26 were sufficient to damp out the torsional oscillations to a 
substantial degree (see Figs. 34-37). It is interesting to note 
that most of the beneficial effect of the feedback is obtained 
with Kp = .13; increasing K p to .26 produces only a marginal 
further improvement. 
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SECTION 5 
CONCLUSIONS 

From the preceding calculations and experiments# the 
following conclusions can be drawn: 

1. Treating stall flutter as a phenomenon induced by 
time-varying aerodynamic damping serves as an 
effective point of departure for analysis and design 
of a stall flutter suppression system. 

2. A relatively simple linear model of the blade and 
servomotor dynamics can serve to give substantial 
quantitative and qualitative guidance as to feedback 
gain selection for a stall flutter suppression system 
based on I.B.C. techniques. 

3. A system centered around the concept of providing 
rate feedback to the blade control motor from blade- 
mounted accelerometers can generate increased damp- 
ing of first torsion oscillations induced by either 
externally imposed or aerodynamic excitation. 

No apparent fundamental obstacle exists to extending 
the control techniques developed herein to larger 
scale systems. 


4 . 
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TABLE 1 

DESCRIPTION OF THE ROTOR BLADE USED 
IN WIND TUNNEL TESTS 


No. of blades 

Radius, w/o steel leaves 

Radius w/ steel leaves 

Chord 

Section 

Lift-Curve Slope 
Drag Coefficient 
Rotational Speeds , Forward 
Flight Cases 
Aerodynamic Center 
Hinge Offset 
Built-in Angle of Twist 



1 

2.03 ft. 

2.311 ft. 

2.0 in. 

NACA 0012 
5.73 
.012 

6.7 Hz (first test) 

6.1 Hz (second test) 
25% chord 

2.0 in. 

8 deg. (linear) 
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FIG. 2 Typical Torsion Loads on Blade Undergoing Stall Flutter 
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FIG. 3 Approximation for Generalized Pitch Damping 
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FIG. 



Side View of Pitch Actuator and Blade Assembly 
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r ; 


FIG. 5 3/4 View of Pitch Actuator and Blade Assembly 



FIG. 6 Individual Blade Control Experimental Rig, Upstream View 
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8 Detail of Blade and Accelerometer Installation, 
"Electronic Spring" Configuration 
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FIG. 9 Detail of Blade Installation, Mechanical Spring 
Configuration 
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FIG. 11 Root Locus Diagram, "Electronic Spring" Static Test 



VoH.k i*p«' 
9>+cV C©** 

C a»<«i 
A«**p 

Qfllifl 

ZI.2? 




«'-*) 



R « A » f 

- .HS I 


Pot 
f 

— I K: - 


r<cib«vL 



?.i 

J \M *\1 ’ 

1 

V/ 1 


5 1 


_l£l_ 

I of 

(milk \o«» " 

roll c?? ) 


Fc«<Ab« ck 
Pot 


0%i V 

(<MW - /-J *•»' 

.0S W ■ .OOlifU 

{ ro^oti^c. ) 


^ccc\e<■u - 
vncter S 


-057 


0.1 ; Jl- C 


»<*= 

= .11 


\(^ r: O.O y O l\ 0.7fc 


0 «« i V i ocvs of H.e f eedljut^ ^aifl Up ft^d tV.e. -VoVa | 


OpS *' l^op Iv b C * rt " H\t Sc * rr e c,S | /\ t " I ^ . / 

(S<?« A» ? pco^u 2 *•«#• a«tw.\i «:(• ts.J 

IG. 12 Control System Block Diagram for Mechanical Spring 
Tests 


















^ roJ/«C J 


Tr.plc *ero oV Ofiyirt 
L omW? ^oVfi q"V * D 


13 Root Locus Diaqrams, Static Tests with Mechanical 
Spring 




K e =.^s 


- 40 - 


VOLTS VOLTS 


PITCH RATE, KP-0. 0 







0 . 


0. 1 


0 . 


T 


T 


2 0. 3 0. 4 

TIME CSEC) 


0 . 5 


0 . 6 0 . 7 


FIG. 16 Closed-Loop, Static "Electronic Spring" Tests 
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FIG. 17 Fourier Transforms of Static "Electronic Spring" 
Responses 
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FIG. 18 Open-Loop "Electronic Spring" Tests, = 6.2 Hz 
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FIG. 19 Closed-Loop "Electronic Spring" Tests, = 6.2 Hz 
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FIG. 20 Fourier Transforms of "Electronic Spring" Responses, 
Q - 6.2 Hz 
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FIG. 21 Open-Loop, Static Mechanical Soring Tests 
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FIG. 24 Fourier Transforms of Static Mechanical Spring 
Responses 
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High Feedback Mechanical Spring Test, u = 6.7 Hz 
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FORWARD FLIGHT PITCH INPUT 
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FIG. 28 Open-Loop, Mechanical Spring Test, H = 6.7 Hz, u = 0.30 
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FIG. 29 Moderate and High Feedback Mechanical Spring Tests, 
fl = 6.7 Hz, u = 0.30 
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FIG. 30 Open-Loop, Mechanical Spring Test, ft = 6.7 Hz 
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FIG. 34 Moderate and High Feedbacks, Mechanical Spring Test, 
ft = 6.1 Hz, u * 0.33 
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APPENDIX 1 

ACCELEROMETER PLACEMENT AND SIGNAL EXTRACTION 

Consider four accelerometers to be rigidly attached to 
the root of a rotor blade, each a distance a from the pitch axis 
(.25 c) as shown below: 



foUt.ow arrows 

When the blade moves about its pitch axis, each acceler- 
ometer will produce a signal purely proportional to a-* if the 
chord of the blade is perpendicular to the axis of rotation 
(i.e. 6=0). However, if 6 * 0, the situation pictured below 
deve lops : 
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Fore accelerometer senses: a 0 - = 

\ coS S \ 

A 0 “ -ft- ^2, ^ — — j <; i*\ 0 — a 0 — _fl ftcoi0ti«0 


Top accelerometer senses: a 0 -*• J1 cT = 


a 0 ♦ _TL 


tX( 


*T. 


cos & — a 0 + _ft2_ cos G 


Sun Signals from Top and Fore accelerometers to yield 2a9; 
purely proportional to 9. 

Exactly the same result could be obtained by summing the signals 
from the Aft and Bottom accelerometers. 
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Note that if the distances of the accelerometers from 
pitching axis are not precisely equal the sum of the 
signals will contain a component proportional both to 
to the difference in the distances. 


the 
two 
9 and 
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Appendix 2 

DEVELOPMENT OF SYSTEM DYNAMIC EQUATIONS 


For the case in which the servomotor shaft and the blade 
have a rigid connection, the motor can be described using the 
following equations of motion: 


il* ♦IR ♦ k E e fc = V (A2.D 

-K t I v i T e c * oe t = 0 (A2.2) 


where 


I = current through motor, amps 

V = voltage input to motor 

6 = angle of motor shaft, radians 

c 

. -2 

L = motor inductance = 1.09 x 10 ohm-sec 

3 

R = motor resistance = 4.93 ohms 
K t = torque sensitivity = .191 n-m/amp 

K_ = volts back EMF = 20 v/1000 rpm = .191 volt-sec 

h 

D = motor viscous friction constant = 3.82 x 10 ^ 
n-m-sec 

J T = sum of motor, tachometer, linkage, and blade 

-5 2 

inertias = 8.35 x 10 kg-m 


The equations of motion can be reexpressed as 
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“ - 


(* L V0 V 
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V 

(A2.3) 
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OJ 

-\< T T t ''-*'Os 


©c 


0 




L J 





To find the transfer function 9 C /I , take the quotient of 9 /V 
and I/V. 




s Vi \ 
-'v o 
MO 

\ V 

o i T o.Oi 
M 0 


v , (A2.4) 

'R 

MO 


V(t t >»0) 

MO 


(A2. 5) 


A(s) = system characteristic equation 


So , 



kr /ft * 
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This, along with the block diagram of Fig. 7, permits the 
derivation of Eq. 1 in the text, so long as the high frequency 
dynamics of the pitch flexure are neglected; if these dynamics 
are to be included, the closed loop transfer function for the 
pitch control system becomes, for K* = .05, K : = .234: 



.T-IS 

(I - 


(A2.6) 


p : =-182.7 + 416. Oj 


p, =-164.2 + 328 . 5 j 


Equation 6 allows the derivation of Eq. 2 in the text. 

To include the effects of the mechanical torsional spring, 
the equations of motion of the blade/servo system become: 


i u* ♦ u ♦ = V 


(A2.7) 




( A2 . 8) 


V + (e* e ‘/O = 0 


(A2.9) 


Equation 9 is the blade/spring equation of motion 
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K = non-rotating torsional spring constant of mechanical 
NR 

spring = 3.53 n-m/rad 

J = inertia of blade, clamps, and counterweights 
B 

= 9.31 x 10“ 5 kg-m 2 

J T = inertia of motor, tachometer, and linkage 
= 8.11 x 10“ 5 kg-m 2 

Note the effects of the 2:1 reduction gear between the servo- 
motor shaft and the blade on Eqs. 8 and 9. 

Equations 7, 8, and 9 can be reexpressed as 



Using matrix manipulation methods similar to those used in the 
previous case, we obtain 
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(A2. 10) 


5000 
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^•0 

l( 

>- 5 40 

'(i- 



‘'UV”) 


MO K„n \'0 0 1 


P 3 = -.057 + 220. lj 

Using this and the feedback gains K = .528 and = .19, 
the closed loop transfer function for the servo block in 
Fig. 12 becomes 


a 

11/ 

( '(wni 

‘ ♦ l] 


v 'l 

f Jz *|)| 

^|W.O / 


s/ fJ 

l" %) 


(A2.ll) 


P 4 = -2.75 + 192 - 9 j 


Including the effects of blade and feedback dynamics as 
shown in Fig. 12 yields the open loop transfer function in the 
following form; 
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_ .1X7 K, t - * 

( '- V f«) ( I * ^.o)' (A2.12) 

where K p includes the gains of both the accelerometer block 
and the integrator block. 

This open-loop transfer function was used to generate 
the root locus diagram in Fig. 13. Loop gain K in that diagram 
is equal to .127 Kp. 

To model the effects of aerodynamic damping on the above 
results, a term D,9 can be introduced into the blade/servo 
equations of motion, Eqs. 7, 8, and 9. The new equations are, 
in Laplace-transformed form. 


(s L *fc»l) I * (“'/nK = ( V) v 


(A2. 13) 


VHJi *( HT K 




»•* 


IQ/ i 




(A2. 14) 


(- '*■) 6c * (Vic/ * ’ l)e = C 


(A2. 15) 


Again using straightforward applications of Cramer's Rule, 
these equations can be used to find 9 c /I: 
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(A2.16) 


As noted in the text, the value of changes with 
azimuth; a "worst-case" value for D^, corresponding to 
$ e ff = -0.1 in the blade equation of motion (Eq. 15), might be 
D a = -3.6 x 10 3 n-m-sec. Inserting this value in Eq. 16 yields 



SOW ‘ ’ l _) 
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p- = + 19 . 0 + 


220 . 9 j 


(A2. 17) 


ror K = . 528, K • = .19, thi.. leads to a closed loop 
transfer for the servomotor of 
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. 7*3/ | 
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p 6 = -*-16 . 9 + 192 . 3 j 


(A2. 18) 


while the open loop transfer function of the entire control 
system becomes 



.IllKfi’ 




(A2. 19) 
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Using Eq. 19, the root locus diagram for the system 

including large negative aerodynamic damping was generated in 

Fig. 13; though this model of the effects of dynamic stall on 

the blade is clearly an approximation (since, in reality, stall 

• •• 

effects have a highly complex dependence on 9, 0, and other 
variables) , the above example shows that the system proposed 
herein can be expected to perf jrm adequately even under very 
adverse conditions at constant negative damping. On a real 
rotor, of course, the blade is exposed to effective negative 
damping only for a small fraction of each revolution. 

Note also that the above analyses have been performed 
under the assumption that the effects of rotor rotation on the 
torsional spring constant K r , R were negligible. In reality 

K R - K NR * J B ^ 

A typical value of ft in these tests was 6.7 Hz, which 

yields 

K r = 3.52 + (9.31 x 10 -5 ) ( (6.7) (6.28) ) 2 = 3.69 

The difference between K R and K^ R was deemed tc in fact be 
negligible for present purposes, especially since the change 
effected due to rotation in the spring natural frequency was 
less than the resolution of the Fourier transform routine used 
in data analysis. 


